(NASA-CB-1 57736) OPTIMIZATION OF MLS 
EECEIVEES FOB MULTIPATH ENVIBONMENTS Annual 
Report {Virginia Dniv.) 62 p BC A04/MF A01 

CSC1 17G 


N78-3208G ) 


Bnclas 


r 


G3/04 31672 


RESEARCH LABORATORIES FOR THE ENGINEERING SCIENCES 



SCHOOL OF ENGINEERING AND 
APPLIED SCIENCE 


40 ^' 

/WN & 



UNIVERSITY OF VIRGINIA 


Charlottesville, Virginia 22901 


Annual Report 

OPTIMIZATION OF MLS RECEIVERS FOR 
MULTIPATH ENVIRONMENTS 

NASA Grant NSG 1128 

Submitted to: 

A. Scientific & Technical Information Facility 
PO Box 8757 

Baltimore/Washington International Airport 
Baltimore, MD 21240 

Submitted by: 

G. A. McAlpine 
J. H. Highfill III 
C.-P. J. Tzeng 
Ghassem Koleyni 


Report No. UVA/528062/EE78/106 
October 1978 



RESEARCH LABORATORIES FOR THE ENGINEERING SCIENCES 

Members of the faculty who teach at the undergraduate and graduate levels and a number of 
professional engineers and scientists whose primary activity is research generate and conduct the 
investigations that make up the school's research program. The School of Engineering and Applied Science 
of the University of Virginia believes that research goes hand in hand with teaching. Early .in the 
development of its graduate training program, the School recognized that men and women engaged in 
research should be as free as possible of the administrative duties involved in sponsored research. In 1959, 
therefore, the Research Laboratories for the Engineering Sciences (RLES) was established and assigned the 
administrative responsibility for such research within the School. 

The director of RLES— himself a faculty member and researcher-maintains familiarity. with the 
support requirements of the research underway. He is aided by an Academic Advisory Committee made up 
of a faculty representative from each academic department of the School. This Committee serves to inform 
RLES of the needs and perspectives of the research program. 

In addition to administrative support, RLES is charged with providing certain technical assistance. 
Because it is not practical for each department to become self-sufficient in all phases of the supporting 
technology essential to present-day research, RLES makes services available through the following support 
groups: Machine Shop, Instrumentation, Facilities Services, Publications {including photographic facilities), 
and Computer Terminal Maintenance. 



Annual Report 


OPTIMIZATION OF MLS RECEIVERS FOR 
MULTIPATH ENVIRONMENTS 

NASA Grant NSG 1128 


Submitted to: 

NASA Scientific & Technical Information Facility 

PO Box 8757 

Baltimore/Washington International Airport 
Baltimore, MD 21240 


Submitted by: 

G. A. McAlpme 
J. H. Highfill III 
C.-P. J. Tzeng 
Ghassem Koleyni 


Department of Electrical Engineering 
Research Laboratories for the Engineering Sciences 
School of Engineering and Applied Science 
UNIVERSITY OF VIRGINIA 
Charlottesville, Virginia 


Report No. UVA/528062/EE78/106 
October 1978 


Copy No. 



TABLE OF CONTENTS 


Chapter 


Page 


1 INTRODUCTION 1 

Background 1 

Objectives 2 

2 STATE- SPACE APPROACH ' 3 

Signal Modeling ^ 4 

3 THE OPTIMAL RECEIVER 6 

Receiver Structure 6 

LOE, and A 10 

4 INTRODUCTION TO THE SUBOPTIMAL RECEIVER 14 

Derivation of the Suboptimal Receiver 14 

5 TABLE LOOKUP TECHNIQUE 19 

6 SIMULATION STUDIES FOR THE SUBOPTIMAL RECEIVER 23 

Simulation Models 23 

Environment and Baseband Receiver Signal 23 

The Optimal Receiver Simulation 25 

The Suboptimal Receiver Simulation 25 

The Antenna Selectivity Function 26 

Threshold Receiver 26 

Simulation Runs and Results 27 


, SBSaa8a PAGE ** " 


11 



TABLE OF CONTENTS (continued) 

Chapter Page 

Crossing Multipath Studies . 27 

RMS Error Studies 40 

Acquisition Scenario 41 

7 CONCLUSION 53 

Crossing Multipath Study 53 

RMS Error Study „ 54 

Acquisition Scenario Study 54 

REFERENCES 55 


111 



LIST OF FIGURES 


Figure Page 

6.1 Crossing Multipath Scenario. Reference Case: High S/N* No 

Interference 29 

6.2 Crossing Multipath Scenario; High S/N, Moderate Interference, 

Low Scalloping Rate 30 

6.3 Crossing Multipath Scenario; High S/N, Moderate Interference, 

High Scalloping Rate 31 

6.4 Crossing Multipath Scenario; High S/N, Heavy Interference, 

High Scalloping Rate 32 

6.5 Crossing Multipath Scenario; Low S/N, Heavy Interference, 

High Scalloping Rate 33 

6.6 Crossing Multipath Scenario; High S/N, Moderate Interference, 

Moderate Scalloping Rate 34 

6.7 Crossing Multipath Scenario; Moderate S/N, Moderate 

Interference, Moderate Scalloping Rate . . . 35 

6.8 Crossing Multipath Scenario; High S/N, Heavy Interference, 

Moderate Scalloping Rate 36 

6.9 Crossing Multipath Scenario; Low S/N, Heavy Interference, 

Moderate Scalloping Rate 37 

6.10 RMS Error Studies; Low S/N, Moderate Interference 42 

6.11 RMS Error Studies; Moderate S/N, Moderate Interference ... 43 

6.12 RMS Error Studies; High S/N, Moderate Interference 44 

6.13 RMS Error Studies; Low S/N, Heavy Interference ....... 45 

6.14 RMS Error Studies; Moderate S/N, Heavy Interference 46 

6.15 RMS Error Studies; High S/N, Heavy Interference 47 


iv 



LIST OF FIGURES (continued) 


Figure Page 

6.16 Acquisition Scenario. Reference Case: Steady-State 

Tracking; 0 se p =1.88° 48 

6.17 Acquisition Scenario; Interference Acquisition when Initial 

Error = 0° and Qsep = 1-88° 49 

6.18 Acquisition Scenario; Interference Acquisition when Initial 

0 R Error = 0.38° and 0 S ep = 1-88° 50 

6.19 Acquisition Scenario; Unsuccessful Interference Acquisition 

when 9 S ep = 1.0° 51 


v 



LIST OF TABLES 


Table 


Page 


1 38 

2 39 


vi 



CHAPTER 1 


INTRODUCTION 

This work is part of the project, "Microwave Landing System " 

(MLS) , undertaken by the Communication Systems Laboratory, Department of 
Electrical Engineering, University of Virginia, under Grant NSG 1128. 
Current work m progress is concerned with the reduced-order receiver 
(suboptimal receiver) analysis m multipath environments. In this chap- 
ter the origin and objective of MLS will be described briefly. Chapter 2 
and Chapter 3 will be the review of signal modeling m MLS, the optimum 
receiver structure, and its performance. Readers are requested to refer 
to the prior reports submitted by the Communication Systems Laboratory 
[1-4] . Chapter 4 will be a summary of the derivation of the suboptimal - 
receiver. Chapter 5 is the description of a computer-oriented technique 
which we used m the simulation study of the suboptimal receiver. 

Chapters 6 and 7 present the results and conclusion obtained from the 
research for the suboptimal receiver. 

Background 

Since man learned how to fly, there has existed a need for a land- 
ing guidance system to aid the pilot during periods of restricted visi- 
bility. The Instrument Landing System (ILS) , which was adopted by the 
International Civil Aviation Organization (ICAO) m 1949, is presently 
the international standard. The limitation to ILS, such as 


1 



2 


susceptability to interference and weather degradation . shortage of fre- 
quency channels, large size of antennas, and the restriction to one nar- 
row approach path has raised the need of a new universal approach and 
landing system. In 1970 the Radio Technical Commission for Aeronautics 
recommended the development of a universal microwave landing system m 
1971. At this time, the United States selected the Time Reference 
Scanning Beam (TRSB) as its choice for the ICAO program. 

In this report MLS will be referred to as the MLS System, i.e. TRSB 
system which has been selected by the United States. 

Obj ectives 

The Microwave Landing System provides an electronic guidance in an 
air terminal area for an approaching aircraft to compute its position m 
space relative to a fixed ground reference. The required information is 
derived by the aircraft's receiver from ground-transmitted microwave sig- 
nals. The goal of the project is to develop an aircraft receiver which 
can give optimal performance m the multipath environments found m air 


terminal areas. 



CHAPTER 2 


STATE-SPACE APPROACH 


State-space approach is the focus of modern fcontrol theory. Sev- 
eral factors influence the development of modern control theory: 

a. The necessity of dealing with a more realistic model of 
the system. 

b. The shift in emphasis towards optimal control and optimal 
system design. 

c. The continuing developments m digital computer 
technology. 

d. The shortcomings of previous approaches. 

State variables consist of a minimum set of variables which are 
essential for completely describing the internal status, i.e. state of 
the system. Conventional input-output equations, or the transfer func- 
tions for linear systems , do not give us any information about the 
internal properties of the system. Optimal control makes it even more 
difficult to avoid dealing with unsatisfactory nonlmeanties , which are 
very difficult to represent m conventional input-output equations. The 
development of modem digital computers makes possible the solution of 
problems which were previously msolvable. Since computers work m the 
time domain, it is more efficient for a computer to directly integrate 
differential equations than to use transform-inverse transformation 
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methods that were usually used m conventional control systems . These 
factors thus justify the use of the state-space approach of modem con- 
trol theory — particularly, as it applies to MLS. 


Signal Modeling 

The whole system is modeled by the state-space approach. The angu- 
lar coordinate [1] to be estimated and other relevant quantities that 
evolve are assembled into an N-dimensional state vector modeled as 
the solution of a suitable linear difference equation evolving m dis- 
crete time, from scan to scan, and excited by a white zero-mean random 
process, | ^ = h -| . 


= (°j 9 , 0s, , 




( 2 . 1 ) 


T • d ( ) 

where ( ) denotes uranspose and ( ) denotes and 


<y = & (k) = direct path signal-to-noise ratio (SNR) 


c ? = e (k) = angular coordinate of own A/C 


o' = cx' (k) = multipath SNR 
R R 


0 = 0 (k) = angular coordinate of reflector specular point 

R R 


(2.2) 

( 2 . 3 ) 

( 2 . 4 ) 

( 2 . 5 ) 


j2 = E (k) = direct-path to multipath phase difference at the 
beginning of the scan 


( 2 . 6 ) 


J3 C l< t t) = 3 (i<y • f- C0-, : 


( 2 . 7 ) 
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where T, = time interval 
k 


(jJ = the scalloping rate 
sc 


( 2 . 8 ) 


The system difference equation can be expressed as 

f(H)x(K)f &Cr<)?(X) 


<1 07 


H. 7 1 f Ky /- / K J j 


(2.9) 


vhere u(k) is the observation, n (k) is receiver noise, 
til 

The 3 component of u, u_^ , can be expressed more specifically as 


j " Jet-pie - (Tj )J -tdp.p [Qh - fa ftj yj cc£ (&+ (jsc. fj ; +■ ncj } 
f ^fj) -f n -'j ] f 


( 2 . 10 ) 


m terms of a discrete-time variable, , local to the scan and, assuming 
the presence of a direct-path component, a single multipath component and 
receiver noise where 


0 ( *) = the transmitting antenna-scanning function 


( 2 . 11 ) 


p [ • ] = rhe transmitting antenna selectivity function 


( 2 . 12 ) 


and 


n , n = independent Gaussian random variables with mean-zero 
<=3 S3 


variance 0.5 


(2.13) 



CHAPTER 3 


THE OPTIMAL RECEIVER 


This chapter contains the summary of the optimal receiver struc- 
ture, operation, and performance. Readers are referred to [2] for 
details. Theory and results m this chapter were used m the next chap- 
ter for the derivation of the suboptimal receiver. 


Receiver Structure 

The objective of the desired MLS angle receiver is to produce an 
estimate of the A/C angular coordinate ,0 , which is minimally affected by 
multipath interference . Recursive state estimation was used in the 
receiver system. If we define 

. } 

y /VJ k ~ 0 l - - - f i = the sequence of observations from 

I ' ' ' ' ' J (3.1) 

some initial time through the present^, 

and 

= estimate of x(k) , given U(k) (3.2) 


then the estimation evolution is described as follows . 



I- ; 



(3.3) 
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where 

7(\<l K~'J r 

estimate of the error in ^ V ^K‘?K"' I J , given Ufh’J- 


(3.4) 

(3.5) 


This is complicated to compute. In our research we assume £ is small m 

£ 

some sense and use the following equation with good results. 

%(m) = 


(3.6) 


where 


estimate of the error m -p.f'A in the neighborhood of zero 

(3.7) 

error, given u(k); 

t y (.K) ~ a gam matrix, depending on y-fK\\-~0 , Q ’><), and statistics of 

eO-4*)> ‘ (3 - 8) 


‘ ~ f Q (■ t Oz 2^ 0 « ; (3.9) 


(Diag. means Diagonal) 

The calculation of g(k|'C) is based on the Locally Optimum 
Estimation (LOE) criterion of Murphy [5] . LOE, as well as the recursive 
estimation Kalman filter, constitutes the following structure of our 


receiver. 
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Fig. 3.1 

As shown m Fig. 3.1, a Locally Optimum Estimator (LOE) used the 
last estimate and extracts all usable information from the new scan data. 
A Kalman filter integrates the output of the LOE with the past to produce 
an optimal estimate, given all data through the present. 

Generally , the Kalman filter uses the following formulas recur- 
sively to achieve the optimum estimate. 


x(ki ,t~ 0 = r { y -0 9“ (V; \ -l 


(3.10) 


p(KIH~>) = T-Cx-0 F(k‘~i itf-i) r ' (y-0 r {•{tf-OG 0 Gt'(h-i) ( 3 . 11 ) 


K'rj - H r - 


^ / ♦ it 
■-1 , i t J f 
* - < ^ I 


(3.12) 


> ' ^ /r J 


S Jr -0 + /<T< J.* /f - J r 1 - H > f .< ' >J " >, 


(3.13) 


fpou Q 





ORIGINAL EASES IS 
03 EOOK QUALKX1 

PiHik) z? P(wik-i) ~ K( yr) HFt K[h*->) 

f 7 H 1 P(h'lk-t) [T-kau !4j r K{*) &lkJ y f h ' 


(3.14) 

(3.15) 




where ^ is a "pre-estimate" of parameter vector using the LOE, e, as 
follows : 


v* (Hi r) - y'LR'.x-O +■ e fJO 


r H5' (‘-I K - >J t r r v 


(3.16) 

(3.17) 


G is as defined m prior work [2] . 

(3.18) 

(3.19) 

where 


% 






< j . . * F 


e f.Vy 




H s masking matrix associated with the choice of V (3.20) 


and, finally, the LOE estimation error is 
y(k) == eiCKj- e ( H) (3.2i) 


Following Murphy's concept [5] m LOE, assuming locally optimum 
estimation, and using LOE pre-estimate, Y , we could simplify the algori- 
thm of the Kalman filter. Equation (3.1J) could be rewritten as [2] 



Jffclvt-i) J- H(t)A(«|^ 
1 ' 0 


(3.22) 


where 



DEK^msi! PAGE IS 
POOR QUALITY 
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M(k) - POUl'-i) H t (T rgHPH 1 ) 


-I 


(3.23) 


th 


and (ti . 1S a vector with i component, as follows: 


a ft?'- i ( ( ^HCy) 

-A..* (ii| % } = j 5 c* ~ ; ' j 

i 0 


if \ / 0 
etherw U.B 


(3.24) 


m which X (',1)1)13 the likelihood ration [3] and 


~ { j\~ -A~ r i t/c tN y 


(3.25) 


(3.14) could also be written as 


PfVjKj 



rt ^ /r F PC^l'^'OC 


N 



•X 


i n 


T 


(3.26) 


LOE,g , and .A 

The concept and development of LOE was expounded by Murphy [5] m 
1968 and summarized m [1] and applied to the A/C angular coordinate 
estimate problem in MLS. 

The LOE first assembled a selected subset of the state vector into 
a parameter vector 7" and then processed the observations u to obtain an 
estimate e of the error m the current*/ - estimate. 


-,*N 


A 



where 


(3.27) 


j£ i. C u > t * J'~ r ('■ I tj l 


VO \N % y 


(3.28) 
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f* t cr) ! 


v - 'i" 


(3.29) 


$ j (y ) ~ U*', I ncj - r.sj — Q , (/> U-CY'i/ 4 j is given by (2.10) 


(3.30) 


where 


y . ^ ^ -j 

-A.(ui^} - l ^ U ^^~ ir 'Jy- £ -, <7 -' 5 


*/ 


(3.31) 


w‘}> '■-fC 


^ . 9 , . r B , ;9 


(3.32) 


and \l'X\ 13 the likelihood ratio, as follows, for J samples/scan: 


a(*!°) = //'A, - T ! 


(3. 32a) 


and I ( ) is the modified Bessel function of the first kind, zeroth 


order. 


Following [2] , we can write 


_A. - 9 '^ 


(3.33) 




(3.34) 


where 


®GTOL T PAqe m 

soon quality 


p ~L ( <? a'* 

u " * ir? rf 




- ■ i A mTRiX 

Dr / ^ 


(3.35) 
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W(u.\ Q.) = 


*j- U; Hi { 4-° »* '-try — I l 
- fu 1 r<r J 7 ' 1 


a j va TiJJ? 


(3. 36) 


m which 


K'i = iL 

dZ 


(3.37) 


M, C*j $ T 


(3.38) 


and 


Hv/£|) = {'v(.ifz}W 7 (i'Jpj f is wtfAsuf 

= D'^ i> (p f W s Cp , U £ f ) | 


(3.39) 


(3.40) 


where 


U L[ (o) d ^ u/j 7 ( u i. ) ! ^ is ionhoof •€■ YYcr^y.- 


( ££* ref. Cv ) 


(3.41) 


s + *?,- 

7 ‘ 


mf OO^Q»lut|3.42) 


th I 

An approximation for the 3 component of w(u|q) in (3.36), above, was 
described m [4] , as follows: 



13 


W j( u l% ) - 


u; 


j \ j - 


S-» J 


( 3 . 43 ) 


°Kj GlhT 





CHAPTER 4 


INTRODUCTION TO THE SUBOPTIMAL RECEIVER 


The optimal receiver generally outperformed the threshold receiver 
at the expense of complexity. Five parameters,^ , Q / and p. , had 

to be acquired before the multipath signal could be tracked m the Kalman 
filter and LOE system. The difficulty m acquiring 0 has. m fact, pre- 
vented the application of the optimal receiver. Consequently, we 
designed a reduced-order receiver, also called the suboptimal receiver, 
which resolved this difficulty while simplifying the receiver structure. 

Derivation of the Suboptimal Receiver 
In the optimal receiver the likelihood ratio, X, involving the 
quantities A- and j? , was used for the LOE system. In the suboptimal 
receiver we were concerned with a likelihood ratio obtained by averaging 
2 out of \ . The parameter vector estimate, V, became a four-dimensional 
vector 

y" - (d g ^ mmmi mm n <4.u 

- J J ' J pg EOOR HgS&BEEe 

Following the derivations in [5], the likelihood ratio, \ , m which |3 is 
averaged out, can be written as 


A ' 


, \ a 

a; = 


T “ 

TT A- 

r* — t ** 


Y", U 


(4.2) 
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>a (X, = j f Aj (V, “j jS S )d ej < 4 - 3 

The structure of the suboptimal receiver is similar to that of the 
optimal receiver, except for the following: 

1. The dimension of the state vector of LOE and Kalman 
filter are less m the suboptimal receiver. 

2. The computation of _A.andjg, which were used m LOE, are 
„ different, since A is a function averaged over 

Formulas for the computing of j\_and 3T were m [5] and are summar- 
ized m the following: 


s>_ \ > 


( S A ? / v* " C’ ^ *j t: ' 

jc. jp ' ' rJ 

ff'Xj u.; , ft'jdG »' 


(4.4) 


-t L 


K ijd 

3r 


/A j 




J ,->» 


where is defined in (3.32a) 


(4.5) 

DKTGIKSD page is 

OE POOR QUALITY 


V 


y Z fj l Co)± 2xa(ft f, ( 9 ) f f ( -r oG P,~ten)~ : 2j 0^4.6) 


If we define 





(4.7) 
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~ Pj ( 1 * 3 


(4.8) 


— fj (&) +" fj ~ (- 6 fiJ 


C | AJ *<>) 


(4.9) 


2J19C&& is) f; (9-) 

\«J w V 


(4 - 10) 


Ej = 


‘T ?A. = 0 

" </ t* « 


('IB;/ i i) 


(4.11) 


Then, _/L.j can written as 


f''' a V< ^ £ 

; P ~ "Fj 

ff- Aj d S‘ 




u &“*&**! 


Jo P Pj 


(4.12) 


So, we can then write 


— P/- \//\ — Pq V/^ 


original 

°$ POOR ( 


(4.13) 


Da- / 


r A . 


-S- Q 
■at-- r J 


Pn -I 


JL. <£ ' 

^ ' j 


: D' 


(4.14) 


and 



17 


VJA ~ C ' 


' Wa J 


• • ) T , Wi * ( 


Vsj ■ - ) 


r 


(4.15) 


and 


(* iiL' cf g ’ 
v-//,j - <? r J 


X* a j «• £j 


(4.16) 


We, 


__ Jo 5 J ; > 4 PJ 

4'>m 


(4.17) 


Also , it follows that 


CP 


- < - A. -A - 1 I Vk ~ r /> 


(4.18) 


= DaHwaDa 1 4- Da biwAS i-V - (& Pf?') T ~r P& H^g P$ (4.19) 


where 


H, 


SA = 


n / ; 

I™ 3 « V. **-* M * r. Vy T -< ^ 




37 I-m 0 .3 (" •> Wf , 

'J 

Nvv «0 — p£«tf (h-WAS: kw'Ati’, 


K V(3 T J 


W'.vA Gj- 


) 


(4.21) 


(4.22) 


(where Diag. means diagonal) and 
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UwAj ~ i'C i - V/ 

- (W 2 I r * - v > 


( 4 . 23 ) 


( 4 . 24 ) 


k y-'A ej - <J v^j vfey ) y vl * y\ 


( 4 . 25 ) 


The functions w_^, w , Hw^ , Hw^ , and Hw^, resulting from the aver- 
aging, are extremely complicated. The computation of these functions was 
done by using numerial approximation on a digital computer. The table 
lookup technique was used m computing Hw^, Hw^, and is discussed m the 
next chapter. 


^QINHIS 3S 

pi] gpOR CSSfclfiS'S' 



CHAPTER 5 


TABLE LOOKUP TECHNIQUE 


In the development of the suboptimal receiver we need to calculate 

the values of the functions Hw , Hw , and Hw , which were defined m 

A B AB 

(4.20), (4.21), and (4.22). Their complex nature is such that they 

cannot be calculated m real time because of the excessive length of the 

required calculation; for example, it requires 16 minutes to calculate a 

set of three functions for a specific set of q and B on the PDP-11 

ah : 

minicomputer and 24 seconds on a CDC 6400 computer. In a complete simu- 
lation run the values of these functions for 6,500 different sets of q 

Ad 

and B_^ are required. The required time spent limits the practical value 
of this approach at present. 

Two approaches were considered to solve this problem. One method 
is to use simple functions which can be computed quickly to approximate 
these complicated functions; the second one is the use of a table lookup 


technique . 


Both approaches need the values of the functions themselves to pro- 
ceed. We first chose some even-spaced values of q^ and B to generate the 
values of those functions and plot them. Here, it was necessary to use 
smooth interpolation to complete the plots; the interpolation was also 

used m the table lookup procedure- Six plots, Hw vs. q , Hw vs. B, 

A A B 

Hw vs. q , Hw vs. B, Hw vs. q , and Hw vs. B, were made. From 
AB A AB BA B 
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these approximate plots, i.e. plots of equally chosen increments, the 
reference values of q and B were then determined more specifically to 

■Ci 

build more accurate plots. The plots were used to find the approximation 
functions, if any, and to decide the reference coordinate intervals for 
the table lookup technique. 

It was seen, from the plots, that these functions were too com- 
plicated and irregular to be approximated by simple functions. At this 
point, approximation functions were investigated for their potential. 
While they were often found to be close to the real functions , they still 
did not satisfactorily reflect the characteristics of the real functions. 

Several things had to be considered in employing the table lookup 
technique : 

1. Would linear interpolation, polynomial interpolation, or 
other kinds of interpolation be used? 

2. Would some modification of the true functions, e.g. 
square root values, logarithmic values, or exponential 
values , be better m interpolation than the true values? 

3. How would the reference points be chosen? 

In this case it is necessary to extend the one-dimensional interpolation 
to the two-dimensional interpolation, since the functions were of two 
variables . 

It was desired to find the simplest interpolation method whose 

deviation was tolerable. Linear interpolation of logarithmic values of 

Kw and I!w and of actual values of 
A B 
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Bab = 






(5.0) 


were shown to provide deviations generally less than ten percent. 

In defining reference points for the interpolation intervals the 
standard procedure demands large intervals m slowly varying ranges and 
small intervals for rapidly varying ranges. The range and reference 
points finally chosen were 


- .7? j? 13 - -77 (5.1) 

, / < $ r A — 10 " (5.2) 

l' 

and 11 intervals for B with reference points 0.01, 0.02, 0.04, 0.06, 

0.01, 0.02, 0.03, 0.05, 0.07, 0.09, 0.95, and 0.99. Also chosen were 24 
intervals for with reference points 0.1, 0.1778, 0.3162, 0.5623, 1, 
1.778, 3-162, 5.623, 10, 17.78, 31.62, 56.23, 100, 177.8, 316.2, 562.3, 
1000, 1778, 3162, 5623, 10 4 , 10 5 , 10 6 , 10 ? , and 10®. A table of 300 
values was then constructed as the first step m the table lookup 
procedure . 

The next step concerned the search problem. Whenever a set of 
and B was obtained, it was necessary to determine m what interval it 
belonged. Two search methods were tried. The binary search was the 
first- Assuming there was no correlation among different q and B 
values, the binary search was conducted by successfully dividing the 
range into two equal parts. Another method, which was termed "the 
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presearch method,” assumes a positive correlation between successive 
values of and B. The current search started with the previous inter- 
val and was then followed by a linear search. Since and B were gener- 
ated randomly, experiments were necessary to determine which method 
proved most efficient m our simulation runs. The results showed that 
the “presearch method” was the fastest* thus, this method was used m the 
table lookup subroutine. 

The final step was to follow the linear interpolation formula 


t 


6r, J 1 :*>,%-,) 5 1 - ax 




f 


<T 


y 2 - 


a* 


( 5 . 3 ) 


* -y fr*, JO ~i l*i.'y J - T frt.-j.,)! *** 


4' 1 


u 

< 1 * 




to obtain the desired values. 



CHAPTER 6 


SIMULATION STUDIES FOR THE SUBOPTIMAL RECEIVER 


Components of the simulation are: 

1. The environment and baseband receiver signal. 

2. The LOE/Kalman filter recursive receiver structure and, 
specifically, both multipath-adaptive and non-adaprive 
variants. 

3. A representation of the Phase III MLS receiver denoted 
the threshold receiver. 

Simulation studies conducted, which focused on the suboptimal 
receiver, included the following: 

A. Crossing multipath interference and comparison for all 
receivers. 

B. RMS error versus 6<sp - & - &R and comparison for all 
receivers. 

C. Acquisition scenario for the suboptimal receiver. 

Simulation models are discussed first and then results presented. 

Simulation Models 

Environment and Baseband Receiver Signal 

Basically, the environmental dynamics are simulated with a state 
model of the form (2.8), without the random excitation, using the state 
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vector, x, (2.1); however, the observations are generated m absolute- 
amplitude form. So, the full model is 


X(k'fi) = ?X(x) , X(o)*X<p 

v( l< ) ~ Kv hlH), O- j y^C^)) 

where 


( 6 . 1 ) 

( 6 . 2 ) 


Sc = the initial state at the start of the simulation 
o 


( 6 . 3 ) 
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( 6 . 4 ) 


C - = rms value of receiver noise at a point in the I-P channel having the 

( 6 . 5 ) 

same signal amplitude as the demodulator output. 


The parameter cr- is assumed known, being a receiver characteristic. 


W> 


a matrix-valued function of its arguments, ;diich compiles the J 
vector v(k) as one with a representative element - -j_(6.6) 
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where 


'■^gXCHNAIi PAGE -> 
p&EJOOR 1 / 


Vj(X ) = trjzMj 


(6-7) 


and u_^ is as given m (2.9). 

Signal data samples are generated only during sampling windows of 
J/2 samples each., located m the TO and FRO scans, respectively, and cen- 
tered where the centroid of the receiver signal pulses are expected. For 
all runs to date 


J = !30 


( 6 . 8 ) 


corresponding to window width of eight degrees m each semiscan. 


The Optimal Receiver Simulation 

The optimal receiver simulation consists basically of the 
following : 


1. Extrapolation of "x 1 to the present, via (3.10). (6.9) 

2. Scan data processor calculation of j\, via (3.33), and?, 

( 6 . 10 ) 

via (3. 34) . 

3. Kalman filter calculation as follows - 


a - ¥('< j >) f Via (3.11) 

(6.11) 

b. Gam matrix, M(k) , via )3.23). 
A 

(6.12) 

c. Y/VjKjr via (3.22) . 

(6.13) 

d. pfkj fj , via (3.26) . 

(6.14) 


The Suboptimal Receiver Simulation 


The suboptimal receiver used the same procedure as the optimal 
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receiver, but the calculations of^A_and $ were different. Calculation of 
-A- follows (4.13), and that of ^ follows (4.19). 

The Antenna Selectivity Function 

The following antenna selectivity function, V[t>) , and its deriva- 
tive ^p/ , £?,)were used in both the optimal and the suboptxmal simulation 
runs ; 
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m which B, the 3 dB’beam width m degrees, was given the value of one 
degree . 


Threshold Receiver 

Performance of the threshold receiver was included m the data for 
comparison. Reference about the threshold receiver can be found m the 
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paper by R. J. Kelly [6], 


Simulation Runs and Results 

The following four parameters are important to the performance of 
an MLS receiver. 


Vw 


( 

ysc. 


Q 




A 


A 


A 


A 


Direct-path signal-to-noise ratio (dB) 

Multipath to direct-path signal amplitude ratio 
Scalloping frequency (H ) 

Q _ 0 ^ , the separation angle between multipath and 
direct-path direction 


(6.17) 


(6.18) 


(6.19) 


( 6 . 20 ) 


The MLS receivers are expected to operate with S/N ratio of 8 dB or 
higher. Values m the range 8 to 20 dB were used m the simulation 
study. 

Another parameter, p , the initial r-f phase difference between 
direct-path and multipath signals, also affected the results. 


Crossing Multipath Studies 

This scenario began with 



0KC3XNSD PAGE $ 
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( 6 . 21 ) 


d p5e p 

cf'f' 



( 6 . 22 ) 


and ran for 100 scans (approximately 7.4 seconds). Runs corresponding to 
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different values of parameters S/N, ^ , and fsc were made. In this scen- 
ario we assumed all runs were mtitialized m the track mode, i.e. all 
estimated variables produced by each receiver were mtitialized to true 
value. Figures 6.1 through 6.9 show the angle estimation errors of the sub- 
opnmal, optimal, and threshold receivers and the composite SNR as func- 
tions of time and separation angle. It should be remembered that the 
separation angle Srsc? si — 2. 7 T [ ^ ^ -L 

j di 

ORIGINAL PAGE IS 
OS EQOR, CJUALITX 

The key parameters S/N, j> (RHO) (BETA) , fsc (FSC) ,' and the rms errors are 
on the bottom of each plot. Figure 6.1 presents time histories of error 
for S/N = 20 dB, |3 = 0, i.e. no multipath interference. Figure 6.2 shows 
the same case, with the multipath signal half as large as the direct-path 
signal. Note that the performance is better for the suboptimal and opti- 
mal receivers m the case of multipath interference. Figures 6.3, 6.4, 
and 6.5 present the time histories for ^ = 180°, fsc = 500 under differ- 
ent SNR and^> . Figures 6.6 through 6.9 present the cases with 
fsc = 51.3 Hz, ^ = -168°, which produce maximum enhancement by the multi- 
path on the TO scan and maximum cancellation on the FRO scan [7] . It is 
clear that the optimal receiver generally performed the best and the sub- 
optimal outperformed the threshold receiver. 

Tables 1 and 2 summarize some interesting features from Figs. 6.2 


through 6.9. 
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Fig. 6.1 

Crossing Multipath Scenario. Reference Case: High S/N; No Interference 
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Fig. 6.2 

Crossing Multipath Scenario; High S/N, 
Moderate Interference, Low Scalloping Rate 
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Crossing Multipath Scenario; High S/N, 
Moderate Interference , High Scalloping Rate 
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Fig. 6.4 

Crossing Multipath Scenario; High S/N, 
Heavy Interference, High Scalloping Rate 
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Fig . 6.5 


Crossing Multipath Scenario; Low S/N, 
Heavy Interference, High Scalloping Rate 
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Tig. 6.6 

Crossing Multipath Scenario; High S/N, 
Moderate Interference, Moderate Scalloping Rate 
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Crossing Multipath Scenario; 
Moderate Interference, Moderate 


Moderate S/N 
Scalloping Rate 
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Fig. 6-8 

Crossing Multipath Scenario, High S/N, 
Heavy Interference, Moderate Scalloping Rate 
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From Table 1, we saw that the performance of the threshold receiver 
is more sensitive to the multipath to direct-path ratio than the subopti- 
mal and optimal receivers. From Table 2 v we noticed that higher SNR gen- 
erally gave better performance and that the SNR affected each receiver to 
about the same extent. 


RMS Error Studies 

These data are sample RMS values calculated for the $ error pro- 
cesses of the three receivers by taking averages over 100 scans of simu- 
lation results run m environments that are stationary, except the phase 
angle|5. The scalloping rate used swept B over a 271 interval during the 
100 scans, thus rendering the averages taken as sample means with respect 
to both noise and . The first ten scans m each run were excluded from 
the averaging to mimmize transient effects in the computation. The 

results are presented as functions of , parameterized by S/N and P, 

\ ' 

as indicated on the various figures. 

The comparative performances elicited by these tests were not 
entirely satisfying m view of the striking contrasts produced by the 
crossing multipath studies. This outcome may be due, in part, to the way 
aborts are processed m the rhreshold receiver simulation (immediate 
reset of the state estimate to the true state on abort) , but the data 
clearly indicates some performance deficiencies in this environment of 
the optimal and suboptimal receivers, as currently structured. Subse- 
quent re-examination of the results of these runs and others indicate the 
optimal and suboptimal receivers, at small separation angles. 


WgWAlj PA<3g to 

"" EQQS QXJAU'S 
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occasionally, at some point in the 100-scan averaging interval, would 
shift to a false-lock null or lose track completely. This is believed to 
be correlated with the sweeping of <3 over a 2 7C interval during the 100- 
scan averaging process. The performance results, given in Figs. 6.10 
through 6.15, appear to be representative of the optimal and suboptimal 
receivers, as presently structured. Elimination of the problems noted 
may require some restructuring of the receiver, possibly reverting ro the 
"nonadaptive ' design when an interference pulse is not both present and 
distinct from the direct-path pulse . This problem is under consideration 
presently. 


iGINAL i'Aw 


Acquisition Scenario £gaa QUAL 1 l v 

There are five parameters ,oi ,9 , , p w , and m the parameter 

vector Y* which is to be estimated m the optimal receiver. If the multi- 


path signal occurs while the receiver is tracking the direct-path signal, 
acquisition of the multipath signal may not occur. The acquisition of 
the parameter ( B , still unmanageable, prompted us to eliminate by averag- 
ing pout of A; however, m the suboptimal receiver V = ,9 , 4* P , Oft, J 

can be acquired. We studied the acquisition scenario for the suboptimal 
receiver and presented the results in Figs. 6.16 through 6.19. 

All figures present the estimation error time histories , respec- 
tively, m ^ , e c<R‘ anc ^ / with the bottom trace on each, showing the 
time history of the ratio Fig. 6.16 m which the receiver was ini- 
tially tracking both the direct path and multipath signals with £}gep 
— . j oy c was used for comparison. Fig. 6.17 m which the multipath 
signal did not occur until the 26th scan presented similar estimation 
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Fig. 6.10 

RMS Error Studies; Low S/N, Moderate Interference 
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Fig. 6.11 

RMS Error Studies; Moderate S/N, Moderate Interference 
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Fig. 6.13 

RMS Error Studies; Lov/ S/N, Heavy Interference 
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Fig. 6.15 

RMS Error Studies; High S/N, Heavy Interference 
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Fig. 6.16 

Acquisition Scenario. 

Reference Case: Steady-State Tracking; 0 se p — 
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Fig. 6.17 

Acquisition Scenario. 

Interference Acquisition when Initial Error = 0° and 6 S ep = 1-88° 
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Fig. 6.18 

Acquisition Scenario. 

Interference Acquisition when Initial 0 r Error = 0.38° and 0sep ~ 
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Acquisition Scenario. 

Unsuccessful Interference Acquisition when 9 se p = 
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after the 26th scan, as m Fig. 6.16. It showed that the suboptimal 
receiver acquired the multipath signal successfully from an initial 
error of 0°. Figure 6.18 shows another successful acquisition for 
©iep - even when the initial error is -0.38°. Fig. 6.19 shows 

an acquisition failure attributable to the reduction of m this run 

to 1.0°. The loss of acquisition capability with diminishing separation 
angle is believed to be related to the steady-state tracking difficulty 
noted m the RMS error studies above. Solution of the acquisition pro- 
blem should accompany the solution of the prior problem. 
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CHAPTER 7 


CONCLUSION 




Previous results had shown that the optimal receiver was generally 
superior to the threshold receiver by a factor of about 20:1; however, 
its complexity was a distinct disadvantage. The order of the state 
vector and the inclusion of m the state vector caused an acquisition 
problem. The objective m developing the suboptimal receiver was to 
reduce the complexity of the algorithm in exchange for an acceptable 
decrease in performance. 

The results obtained m rhis study can be summarized m three 
sections . 


Crossing Multipath Study 

The integrated LOE/Kalman filter suboptimal receiver algorithm 
tested m simulation was generally superior to the threshold receiver 
but, as expected, inferior to the optimal receiver. The reduction of the 
order of state vector and parameter vector simplified the structure of 
the receiver itself. The employment of the table lookup technique did 
speed up the computation required in the suboptimal receiver; however, 
the length of the computation time still limited the capability of the 
suboptimal receiver, because we used a general-purpose minicomputer. By 
adequate usage of special-purpose microprocessors or good approximation 
functions, this problem may be eliminated. 
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RMS Error Study 

Difficulties m steady-state tracking were noted that appear to be 
analogous to the difficulties m resolving two closely spaced radar tar- 
gets with a finite aperture antenna. "To resolve, or not" is a question 
that needs to be answered (dynamically, as a function of ) m the MLS 

receiver problem, and an answer m the negative cannot be immediately 
ruled out as second best in this case- The problem is being studied. 


Acquisition Scenario Study 

The difficulty m acquiring p was solved by averaging p out of the 
parameter vector m the suboptimal receiver. The simulation study 
revealed both success and failure m acquisition. A more thorough study 
of the acquisition problem is m progress to establish that the subopti- 
mal receiver can acquire when it can track. 
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